Wigner crystals of ions as quantum hard drives 
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Atomic systems in regular lattices are intriguing systems for implementing ideas in quantum 
simulation and information processing. Focusing on laser cooled ions forming Wigner crystals in 
Penning traps, we find a robust and simple approach to engineering non-trivial 2-body interactions 
sufficient for universal quantum computation. We then consider extensions of our approach to the 
fast generation of large cluster states, and a non-local architecture using an asymmetric entanglement 
generation procedure between a Penning trap system and well-established linear Paul trap designs. 
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Quantum information processing using trapped ions 
has been the focus of theoretical [l| and experimen- 
tal @, H, fl S @1 efforts over the past decade. The co- 
herence times of ions can exceed seconds, while manip- 
ulation and entanglement time scales can be as fast as 
tens of microseconds. So far, approaches to scaling these 
systems to many ions have met with significant issues, 
both in linear Paul trap systems, where increasing num- 
bers of ions leads to control difficulties, and in proposed 
more complex trap arrays, where "shuttling" of quan- 
tum information using gate electrodes would allow for a 
scalable architecture [7|. A possible solution is to sepa- 
rate the processing elements (processor qubits) from the 
memory [8(. 

One natural system to consider as a quantum memory 
is a Wigner crystal of ions in a Penning trap [9]. Such 
crystals can be robustly formed [l(| , and are dynamically 
stable, with tens of thousands of ions in a given trap. In 
addition, the strength of the Coulomb interaction leads 
to large separations between individual ions, making in- 
dividual addressing of ions in such lattices a distinct pos- 
sibility, in contrast to present control in neutral atom and 
polar molecule lattices 
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In this Letter we develop an approach to quantum 
memory and entanglement generation that takes full ad- 
vantage of the advances in ion trap technology for build- 
ing large Wigner crystals of ions in Penning traps. Us- 
ing a modulated-carrier "push" gate adapted from linear 
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ion trap quantum computing schemes [13l . IL 
we find a fast but adiabatic method for building small 
clusters of entanglement which is insensitive to thermal 
phonons in 2D and 3D Wigner crystals. We take ad- 
vantage of some of the unique features of Penning traps, 
such as rotation of the crystal, to provide simplifications 
in the necessary hardware to implement these ideas in 
2D Wigner crystals. We further show that such a quan- 
tum memory device can also be used directly for cluster 
state quantum computation. Our approach follows recent 
work [17| on performing quantum gates in 2D Wigner 
crystals. Finally, non-deterministic entanglement gen- 
eration between distant ions suggests a processor (linear 



FIG. 1: Two-qubit gate via intensity modulation of a laser 
addressing a pair of ions in a crystal rotating at frequency ui r . 



Paul trap) and memory (2D Wigner crystal) architecture 
based upon a quantum register approach 0, [l8| , where 
the low photon collection efficiency from ions in the mem- 
ory is offset by an asymmetric entanglement generation 
scheme u sing a weak cavity coupled to ions in the pro- 
cessor 0, [ij] . 

We start by considering a Wigner crystal of ions, ro- 
tating in a Penning trap [9j with harmonic confinement 
with frequencies uj xv (in the lateral directions) and oj z (in 
the vertical direction). With characteristic ion spacings 
d ~ 10 fim, tightly focused lasers allow for individual 
addressing of ions (see Fig. []}. Laser cooling can re- 
duce the temperature ss 1 mK, yielding on the order of 
10 2 — 10 3 phonons in the softest (lateral) modes. By us- 
ing long-lived, metastable states of the ions as a quantum 
memory, we may neglect memory errors in our discussion. 
A tightly-focused laser allows for nearest-neighbor phase 
gates and for single ion operations. Large-scale computa- 
tion may be considered using either nearest-neighbor cou- 
plings or via a variety of quantum communication tech- 
niques developed for quantum repeater protocols. When 
used in conjunction with the deterministic phase gate de- 
veloped below and local single ion operations (performed, 
e.g., via Raman transitions), this will suffice for perform- 
ing large scale quantum algorithms [l8| by using the re- 
mote CNOT gate [H|. 

Modulated-carrier gate — A spatially inhomogeneous 
laser detuned from the appropriate transitions between 
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internal (qubit) states of an ion (a two-level system with 
Pauli matrices a* ,y ' z ) produces a ponderomotive force 
fi due to the gradient in its intensity. Using an appro- 
priate combination of polarizations and frequencies, in 
analogy with alkali atoms [2(| , the sign of the force be- 
comes dependent upon the internal state of the ion, with 
the associated perturbation to the system: 
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where s£( is the displacement of ion i away from its equi- 
librium position. The latter can take place either al ong 
the separation between two individual ion microtraps [l4j 
or perpendicularly to the plane of an ion Wigner crystal 
in a Penning trap [13] • In both schemes, as the ion dis- 
placements are coupled (via phonons), such a push leads 
to an effective erf a? interaction. Adiabaticity is required 
for vibrational excitations to be absent after the gate. 
This bounds the clock speed to be lower than the fre- 
quency of trapping in the push direction: tight traps are 
needed for fast, temperature-insensitive operation. 

We now introduce a simple variant of the fast-kick 
"push" gate which allows us to use even the soft (lat- 
eral) modes when their temperature is extremely high. 
Our variant uses slow modulation of a fast, oscillating 
state-dependent force. The oscillation averages any ion 
motion to zero over the course of the gate, while the 
in-phase oscillation of nearby ions leads to a non-trivial 
phase evolution and the desired quantum gate between 
ions in the crystal. In addition, as our gate allows for 
non-trivial oscillation of ion positions in all three spatial 
dimensions (versus only in the tightly confined direction 
for the vertical gate), it can work using a single laser 
beam and in three dimensional crystals. 

It is instructive to recall the general descri ptio n of 
"push" phase gates when in a complex crystal [l5|, [17 1 . 
We start by rewriting the Hamiltonian of N interacting 
ions to second order in displacement from the equilib- 
rium positions H = ^2 K Hi^k^k^k, using normal-mode 
coordinates indexed by K = {k, A} (the wavevector and 
polarization), f 4 = J2k M iK e K (u K j \[2)(a K + a K ). The The pulse-shape independent form factor is 
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which are exact to second order. 

We now seek an approach which still maintains no net 
change in displacement and no dependence of the overall 
phase on phonon state, but can operate on time scales 
on the order of ujk- We add a sinusoidal variation to 
the force \f[t) — ► cos{vt)f{t)\. The carrier frequencies v 
must be fast with respect to lok', qualitatively, this av- 
erages out any net displacement. If the modulation f(t) 
is slow as compared to v (but with no restriction with 
respect to <jJk)i we can perform a similar adiabatic elim- 
ination as above, and get a gate with the same desirable 
properties that can operate non-trivially on arbitrarily 
"soft" phonon modes at very high tempeatures. 

For adiabatic elimination with respect to v, we choose 
the ansatz — [3+e wt + /3_e~ wt (subscripts omitted for 
clarity). Setting f3 = yields (3± = af{t)/[2y/2h(u±v)]. 
We find the displacement of a normal mode induced by 
the gate is proportional to the force applied, and can be 
made zero independent of initial phonon state by starting 
and ending with zero force. This eliminates any poten- 
tial error due to entanglement between phonons and the 
internal states of the ions. 

We now examine the two-body phase induced in this 
new scenario. The differential equation for phase is now: 
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Averaging the quickly varying component lets us replace 
cos 2 (i4) with 1/2. Returning the mode index, K, we find 
the overall phase accumulated, J2k ^(T), for a gate oc- 
curring over a time to T does not depend on the phonon 
initial state. However, the internal states of the ions are 
affected by the unitary exp(— iYlij ^ij^i^j) where the 
two-body phases are given by 
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uk = \/h/mujK are the oscillator ground state lengths; 
the matrix M is orthogonal {M*M = MM 1 = 1). The 
perturbation V can now be written as 
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V = V a K f K (t)(tf K + a K )/V2, 
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where fx{t) is the state-dependent force on normal mode 
K defined via the transformaton M and Eqn. [TJ 

The problem factorizes into 3iV independent, driven 
oscillators. For scenarios with lim t ^± 00 /(t) = 0, the 
oscillator evolution is given by the unitary transform 
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(the polarization vectors eV only depend on A). 

Expanding in inverse powers of the large carrier fre- 
quency v, we note that the first term, 0(v^ 2 ), is pro- 
portional to J2 k Mik v \Mjk,\ = (due to the orthog- 
onality of the matrix M). The first non-zero term is 
0(z^~ 4 ). Compared to adiabatic gates, this modulated- 
carrier push gate is inverted in sign and reduced in phase 
by a factor (u)/v) A /2. For lateral gates, u> ~ aj xy is a 
characteristic confinement energy for a single ion in the 
crystal, and for vertical gates, ui ~ u> z . 
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Performance — We performed numerical simulations 
of the modulated-carrier gate's performance for finite-size 
2D and 3D Wigner crystals (N = 147 shown in Fig. ^ 
to compare to the equivalent adiabatic gate and the pro- 
posed vertical gate of Ref.[l3- Our simulations minimized 
the classical energy of the ions in a Penning trap to deter- 
mine the equilibrium configuration. Then, expanding to 
second order in displacements from equilibrium, the nor- 
mal mode coordinates were found. Coriolis forces were 
neglected; their inclusion does not qualitatively change 
our results. A gate between the two centermost ions was 
simulated by computing the displacement (3k and phase 
4>k for each phonon mode. The fidelity was calculated 
by considering the overlap of the final state with the de- 
sired state, traced over the phonon degrees of freedom 
and minimized over all possible initial states. 

We find that for the same physical parameters, the ra- 
tio of forces (i.e., laser power) required for achieving a it 
phase for the vertical gate and the fast carrier gate goes 
as (uj z /v) 2 , consistent with the ratio between adiabatic 
and fast carrier gates derived above. Thus, the fast car- 
rier gate requires substantially less laser power for the 
same conditions with negligible reduction in fidelity. Al- 
ternatively, the gate time could be reduced, enhancing 
the overall performance of quantum information proto- 



cols. For specificity, setting uj xy — 200 kHz, 
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MHz, and a gate time r = 5 fis, we find v — 2.2 MHz 
provides 1 — F < 10 -5 with negligible heating. Even 
smaller errors are found in simulations of the 3D crystal 
under the same approximations. 

A practical limitation occurs due to the spontaneous 
emission induced by the off-resonant laser interactions. 
Tight focusing increases the force for the same laser 
power; thus, using a pair of adjacent, narrow-waist (< 
2 /zm) laser beams reduces spontaneous emission and 
power requirements. For specificity, using a transition 
with spontaneous emission of 7 = 20 MHz and lasers 
with peak Rabi frequency of 100 GHz detuned 100 THz 
from the atomic transition, a laser power of ~ 3 mW per 
beam is required for our gate, with an induced error of 
< 0.1% per gate. 

Quantum cluster state generation — We now consider 
an approach that takes advantage of the Coulomb inter- 
actions in the lattice to create and dynamically extend a 
cluster state for universal measurement-based computa- 
tion. Specifically, the goal is to obtain a weighted-graph 
state exp(i^^- afa^'&ij /2) | + , ...,+) , where in the ideal 
case dij equals tt between nearest neighbors on a square 
lattice, and zero otherwise. On a triangular lattice like 
the one available in many-ion Penning traps, this can be 
achieved if dij is made to vanish along one side of each 
lattice cell, and to be tt on the other two. The idea is to 
obtain this via a global tt/2 qubit rotation followed by 
a push gate acting on all three cell vertices at the same 
time, possibly with a laser swept at constant velocity 
through the cell itself, to take advantage of the uniform 
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FIG. 2: (top) Phonon spectrum for vertical and lateral 
phonons for a trap with N = 147 ions and u) z /u) X y = 50. 
(bottom) Fidelity versus temperature for the vertical-phonon- 
mode gate of Ref. [17J (red-dashed line) and for the fast car- 
rier gate with v = llu> xy between the center-most pair of 
ions. Left: fast carrier gate's and vertical gate's forces on one 
of the two ions over the gate time; both gates operate in a 
time r ~ l/uj xy . Anharmonic corrections to the fidelity are 
not included here. For this choice of parameters, the vertical 
gate [TtJ requires 20 times the force (and laser power) of our 
fast carrier gate to achieve the same final n two-body phase. 
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FIG. 3: (color online) a) A laser is swept adiabatically from 
left to right, leading to a weighted graph state with different 
phases for solid and dashed lines, b) Laser displacements in 
the crystal plane needed to obtain a constant sweep velocity 
at given distances £/R = 0.2,0.4,0.6,0.8 (from blue to red) 
from the center, in the rotating crystal's frame. 



circular motion of the lattice. 

We start by considering a focused laser beam of waist 
a (in units of the lattice length d) adiabatically swept at 
constant velocity v through the Wigner crystal, along a 
direction parallel to one of the lattice vectors (Fig. [3]), 
at half the height of a triangular cell. The effect of 
this sweep is, apart from a global single-qubit rota- 
tion, to generate a weighted-graph phase, where in- 
takes value e9(ui) on the cell side that is parallel to the 
sweep direction, and 6{ui) on the other two sides, with 



e -3/(8a )( H _ Sa 2 )/(a 2 + 8), while 
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where a = y h/(mui), firj is the peak Rabi frequency 
corresponding to the center of the laser beam, A is its 
detuning from the ion's internal transition, and q is the 



4 



a) Memory 
(Penning 
trap) 



Processor 
(Paul trap) 



c 




Interconnect 
(beam splitter 
^ + detectors) 




FIG. 4: a) Schematic of a quantum processor (such as a linear 
Paul trap) coupled via a high finesse cavity to a photodetector 
system to allow interference with photons from an ion Wigner 
crystal in a Penning trap nearby, b) Ion level structure 



electron charge. Using the fast carrier modulation de- 
scribed above, the semiclassical calculation of Eq. ([7|) is 
no longer valid, but the discussion of Eq. ^ shows that 
the resulting phase is simply —6{y)/2. A cluster state 
is then obtained by making e small via an appropriate 
choice of the laser waist (numerically, a < 0.2), while 
tuning 6(y)/2 to n by adjusting the other experimental 
parameters such as laser power. 

Care needs to be taken to ensure a sweep having a 
given distance £ from the trap center and velocity v in 
the rotating crystal's frame. To this end we apply to the 
laser, initially focused at a distance R from the center, a 
displacement 5d^(t) = {6x^(t),6y^(t)} of the form 



Sx s (t) = [R-Z/cos(GJ r t)]e(x- Ml) , 
<%(*) = £[wr*tan(x)/x - tan(w r i)]9(x • 
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(see Fig. [3]), where \ = arccos(£/i2). Cluster state gen- 
eration as presented uses a single laser beam, resulting 
in substantially higher laser power requirements than the 
two qubit gate with two beams as described above. In 
particular, for errors per gate < 0.1%, a detuning of 200 
THz and peak Rabi frequency of 4 THz (corresponding 
to 5 W of laser power) would be required. More complex 
laser motion could reduce these requirements. 

Asymmetric entanglement generation — We conclude 
with a brief discussion on the implementation of circuit- 
based computation using entanglement generation and 
remote CNOTs. We use a quantum processor unit (such 
as a linear Paul trap) separated from the quantum mem- 
ory unit (our Wigner-crystal-based quantum hard drive - 
see Fig. 0} , characterized by photon collection efficiencies 
77 and rf respectively. Without loss of generality, we will 
assume 77' > 77, as it can be achieved via coupling with 
high finesse cavities 0, [3] . Our two-click asymmetric en- 
tanglement generation procedure starts with an equally 
weighted superposition I +) = |0) + |1). An optical 7r pulse 
produces photons via spontaneous emission at a rate 7 
from the |1) <-> \E) transition. Then the photons are in- 
terfered on a beam splitter. Without assuming photon 
number resolving detectors, the state after one "click" 
becomes 



(10) 



with \ip±) = \J r\ I T) 1 (|01) ± 1 10) ) . To symmetrize the en- 
tangled state and simultaneously remove the |11) compo- 
nent, a 7T pulse between the metastable states (|1) <-> |0}) 
followed by repetition of the above protocol results in a 
pure state |$+) = (|01) ± |10))/\/2. The overall proce- 
dure succeeds with probability 7777', indicating that the 
time required is (77777') _1 . A standard one-click scheme 
[22] ] with excitation probability p takes a time (77773) 
and succeeds with error rate 0(p), i.e., the higher fidelity 
a pair one wishes to generate, the longer it takes. By 
contrast, in our scheme the fidelity can be high without 
a further increase in generation time. 

Thus, for large-scale computation, a central proces- 
sor unit with high collection efficiency allows for high- 
fidelity gates between elements of the "hard drive" mem- 
ory on a timescale 2/Trjrj' (see Ref. |JLS| for further im- 
provements). For concreteness, we take a radiative decay 
rate of T = (2-7r)10 MHz, 77 = 10~ 3 (confocal approach 
with low numerical aperture lens), and desired infidelity 
1 — F < 10~ 4 . Entanglement generation between two 
such ions would take a time ~ 10 ms or longer; in con- 
trast, for 77' ~ 0.1, using the intermediate quantum pro- 
cessor leads to entanglement generation between proces- 
sor and both ions in a time of order 100 /is, comparable 
to the phase gate operation times already discussed. 

This complements the quantum hard drive architec- 
ture described above, providing a comprehensive toolbox 
for universal quantum computation with ion crystals in 
Penning traps that relies on existing technologies under 
available experimental conditions. 
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